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Abstract. We construct a Hamiltonian formulation of quasi-local general relativity using 
an extended phase space that includes boundary coordinates as configuration variables. This 
allows us to use Hamiltonian methods to derive an expression for the energy of a non- isolated 
region of space-time that interacts with its neighbourhood. This expression is found to be 
very similar to the Brown- York quasi-local energy that was originally derived by Hamilton- 
Jacobi methods. We examine the connection between the two formalisms and find that 
when the boundary conditions for the two are harmonized, the resulting quasi-local energies 
are identical. 
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1. Introduction 

General relativity is a fully covariant theory of gravity and as such does not privilege any 
particular flow of time. This is one of its great virtues but it is also a problem if one wishes to 
study it using traditional methods that manifestly depend on a notion of time. For example, 
a flow of time must be defined before one can do a standard Hamiltonian phase space analysis 
of a field theory. Thus, if we want to apply such methods to Einstein's gravity, space-time 
must be artificially broken into space and time. We do this by foliating space-time into a set 
of space-like three-surfaces ("instants of time") S t , and defining a "time-flow" vector field 
T a that evolves these surfaces into each other. 

Having done this, one can reformulate general relativity in terms of phase space, 
symplectic structures, and Hamiltonians. In the most common approach, the space-like 
three-metric hah is taken as the configuration variable while its conjugate momentum, P ab 
is closely related to the extrinsic curvature of the three-surface in four-dimensional space- 
time (see for example PUEJIS])- Traditionally, space-times studied in this way were either 
boundary-free or taken to have a boundary at spatial or null infinity (with appropriate 
conditions imposed so that the space-time is asymptotically flat). More recently, people 
have become interested in studying general relativity over finite regions of space-time, which 
are often referred to as quasi-local regions. Then, boundaries and boundary conditions at 
locations other than infinity must be considered. In the assorted boundary conditions 

that give rise to a phase space with a well-defined symplectic structure were studied. In 
particular, it is well known that this can be done if one fixes the intrinsic metric 7 a b of the 
boundary to be constant and not affected by variations. That is, one fixes the intrinsic 
geometry at each point of the boundary manifold. 

However, some aspects of this fix are not entirely satisfactory. For example, if the fixed 
boundary metric is not axially symmetric, then a rotation along a vector field <p which is 
tangent to the boundary at the boundary is not an allowed variation in the phase space 
under consideration, because <^7 & = C- <fiab 7^ and so the variation does not preserve 
the boundary condition. Thus, there is no Hamiltonian generating this motion in phase 
space and consequently no notion of the angular momentum associated to such a rotation. 
Similarly, time translations are not permitted unless the boundary metric is invariant in time, 
i.e. Sr^jab = £-Tlab = 0. These evolutions cannot be generated by Hamiltonian functionals 
unless <p a and T a are Killing vector fields of 7^. 

A little thought shows that this is really not so surprising. In both of the examples 
discussed above, the "conserved" quantity corresponding to the listed translation is not 
conserved. Hence, we should not expect to be able to obtain it from a standard Hamiltonian 
treatment, which by its very nature applies to situations where Hamiltonians are conserved. 
For example, if a Hamiltonian exists that generates time translations, then Hamilton's 
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equations are 



5H T = Q(5 T , 5) 



(1) 



where the symplectic structure f2 is antisymmetric in the two variations. Hence, it follows 
immediately that StHt = 0. Similarly, the Hamiltonian generating a particular rotation is 
conserved under that same rotation. 

Despite the arguments given above, Brown and York have introduced notions of 
the energy and angular momentum associated with a bounded region of space-time (if 
the intrinsic boundary metric is fixed) by performing a Hamilton-Jacobi analysis [Hj. 
Furthermore, the expressions they have obtained have proved useful in a very large number 
of applications (as a representative sample see |Sj). It would be very surprising if these 
expressions could not be derived from a phase space treatment of general relativity over a 
manifold with boundary. 

The purpose of this paper is to show that it is indeed possible to obtain these expressions 
for energy and angular momentum using a careful phase space analysis of general relativity 
on manifolds with boundaries. Here, we begin with the standard phase space formulation 
of general relativity in terms of ADM variables (see, for example, jH] for details). The key 
idea is then to import extended phase space techniques from classical mechanics which are 
designed to deal with situations where "conserved" quantities are not conserved. The basic 
idea is to enlarge the phase space under consideration by including quantities such as the 
time coordinate as configuration variables (see ^Hj for a standard reference). It is then 
possible to construct a conserved Hamiltonian generating time translations in the extended 
phase space. However, the energy of the system is no longer the value of this Hamiltonian, 
but is instead the value of the momentum which is canonically conjugate to time. 

In addition to extending the phase space however, it is also necessary to weaken the 
boundary conditions so that time translations and rotations are allowed variations. For 
example, in the Brown- York treatment the boundary three-metric 7 a b is completely fixed on 
a three-manifold B that is a boundary of the space-time manifold M. Therefore, 



for all variations. That is, the metric is fixed with respect to the manifold. Under our looser 
treatment, we allow the variations to act as diffeomorphisms that move our boundary fields 
around B. Thus, 5 acts on the "fixed" fields as an infinitesimal diffeomorphism and so 



for some vector field 5X a . Picking 5X a appropriately, we obtain translations and rotations 
as variations. 

The generalization of the Hamiltonian formulation for a finite region of space-time is, 
to our knowledge, presented here for the first time. However, similar ideas have previously 



Slab 







(2) 



0~lab = £ SXlab-, 



(3) 
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been discussed at spatial infinity. Notably, in their early work on boundary terms for the 
gravitational Hamiltonian, Regge and Teitelboim 2,\ implemented a similar programme for 
the boundary at spatial infinity in asymptotically flat space-times. Specifically they included 
variables that located the asymptotic position of the boundary as well as their conjugate 
momenta as canonical variables and showed that the resulting Hamiltonian was covariant 
with respect to the asymptotic Poincare group. Kuchaf has also studied parameterizations 
at infinity. In his considerations of Hamiltonian formulations for spherically symmetric black 
holes [Zj, he included the Killing time as a canonical coordinate and found that the mass 
was conjugate to its radial rate of change. Finally, Kijowski [5| has also considered general 
relativity in a manifold with boundary. He begins with a novel and non-standard approach 
to symplectic geometry and the Hamiltonian formulation JT] . Despite this different starting 
point, he obtains energy expressions for the boundary which are analogous to the Brown- 
York. However, he only explicitly considers the case where the evolution is a symmetry of 
the boundary. Thus, our work may be considered a generalization to arbitrary boundary 
geometries. It is likely that our method could also be incorporated into this alternative 
formalism. 

The algebra for the quasi-local gravitational case is quite formidable, and so as an 
introduction, we begin in section |21 with a simple example of the extended phase space 
formalism, namely a time dependent harmonic oscillator. With this experience in hand, 
in section El we apply similar techniques to the problem of interest: general relativity in a 
space-time manifold with a boundary. In particular, we obtain expressions for the energy and 
angular momentum associated to the region of space-time, which may be calculated using 
just the values of the fields at that boundary. Section 0] provides a comparison between our 
results and those of Brown and York. We end with a discussion of the results and possible 
extensions and applications of this work. Several key technical results are collected in the 
appendices. 

2. An Introductory Example 

In this section, we shall consider a simple example which will demonstrate how a phase space 
can be extended to allow the description of systems in which the energy is not constant. By 
starting with a time-dependent harmonic oscillator we will capture many of the central ideas 
of the construction without the extra complications that arise in gravity. The key idea will 
be to extend the phase space by including the time coordinate i as a configuration variable 
as well as its conjugate momentum p t . Then if we know how the variables evolve in time, we 
can manipulate the symplectic structure so as to find a Hamiltonian function which generates 
that evolution on-shell. Equivalently, we solve the Hamiltonian equations of motion for the 
given evolution to find the corresponding Hamiltonian. 
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With the standard phase space treatment, the energy of the system is the on-shell 
value of this generator of time translations (which is, of course, constant). In the extended 
phase space treatment, the energy associated with a time translation is equal to the negative 
of the value of the momentum conjugate to t. A slight complication is that, in general, 
the Hamiltonian is not unique and as a result of this ambiguity the energy conjugate to 
the evolution is not unique either. Indeed a significant degree of freedom remains in its 
definition. We illustrate and elaborate on these issues in the following example. 



2.1. Simple harmonic oscillator 

Consider the harmonic oscillator with mass m and spring constant k. The canonical 
phase space of this system is parameterized by the position coordinate q and its conjugate 
momentum p. The symplectic structure is simply given by 

n(5 u 5 2 ) = (5 iq )(5 2 p) - (5 2 q)(5 lP ) . (4) 

In this and all future expressions, one should keep in mind that the 5s are vectors 
in the phase space of all possible configurations of the system. Their conventional 
interpretation as infinitesimal variations of system configurations arises from considering the 
one-parameter families of phase space diffeomorphisms that are generated by such vector 
fields. "Infinitesimal" changes of those parameters generate what we intuitively think of as 
infinitesimal variations of the system. 

Then, as discussed in the introduction, a phase space vector field 8 t (equivalently an 
evolution of the system) is said to be Hamiltonian if there exists a function H t such that 

n(5 t ,5) = 5H t , (5) 

for all variations 5. Conversely, given a Hamiltonian function H t , we can find out how it 
evolves the system by solving the above equation for 5 t . Thus, there is a mapping between 
Hamiltonian evolutions and Hamiltonian functions. 

As an example, for the evolution 5 t which gives rise to the usual equations of motion: 
dq p dp 

— = — and — = -kq, 6 
at m at 

we can show that 

Q(S t , 5) = 5p - (-kq)5q = SH t , (7) 

where 

and C is a free constant. Thus, the time evolution is generated by H t , which the reader 
will immediately identify as the classical energy of the system (up to a constant). Note 
too that the equations of motion confirm that 5 t H t = (which we knew already by the 
skew-symmetry of the symplectic structure). 
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2.2. Time- dependent simple harmonic oscillator 

Things get more complicated if k is not a constant, but instead varies in time. In this case, 
the energy of the oscillator will not be constant, but instead may vary as a changing k adds 
energy to or removes it from the system (more properly, the external agency setting k can do 
net work on the system). Thus, if we have no knowledge of how k is being fixed, we cannot 
do a standard Hamiltonian analysis of the system - such calculations require a closed system 
which does not exchange energy with its surroundings. However, we can use the extended 
phase space formalism to partially compensate for our ignorance of the external system that 
sets k. All we require is that k(t) be a fixed function of t. 

To allow for a time dependent spring constant k and consequently a time dependent 
energy, it is necessary to extend the phase space to include t and its conjugate momentum 
Pt- The symplectic structure is then given by 

n{8 u 5 2 ) = {S iq )(S 2 p) - (5 2 q){5 lP ) + {5 1 t){5 2 p t ) - (S 2 t)(8 lPt ). (9) 

Furthermore, we would like to allow general evolutions, rather than restricting to (d/dt). 
Thus, we shall study evolution generated by 

A = A |, (10) 

where A G is a free parameter which is required to be strictly positive. Our task is then to 
find a Hamiltonian \K A which generates the following evolution: 

p 

S A q = A G — , 5 A p = -X Q kq, and 8 A t = X . (11) 
m 

Note that the evolution of pt under 5 A is not determined a priori. This in turn leads to a 
freedom in the form of the Hamiltonian K A . 

We proceed by evaluating Q(5 A ,5). Making use of (fTTj) . as well as the fact that k(t) is 
a fixed function of t so that 

5k = k 5t 

we obtain 

n(*A, S) = A 5 (vt + + tpj - {^p- + s AP ) st 

= 8{\ K t ) - K t 5\ - + 5 AP )j 8t , (12) 

where we have defined 

% We will follow the convention of ^U] an d use K to denote the Hamiltonian in extended phase space. 



Canonical Phase Space Formulation of Quasi-local General Relativity 



7 



Now, our evolution 8\ will be a Hamiltonian vector field in the phase space if and only 
if £1(8 \, 5) is an exact variation. This will only be true if the 5X and 5t terms vanish. 
Therefore, we must be at a point in phase space where 

Furthermore, p t must satisfy the equation of motion: 

OAPt = — , (15) 

which guarantees that the constraint (fPfj) is preserved under A-evolution. Then, on this 
constraint surface, the evolution is Hamiltonian and generated by 

K A = X K t . (16) 

The reader will immediately realize that this function vanishes on-shell. However, this does 
not mean that the energy of the system will vanish. In the extended phase space, the energy 
is given by the negative of the value of the momentum canonically conjugate to the time. 
Thus we obtain 

E t :=- Pt = ^+ k -l. (17) 

This is immediately recognized as the usual energy associated to a harmonic oscillator, 
although it will not necessarily be constant due to the time dependence of k(t). 

The Hamiltonian K\ given in ()16|) generates the desired evolutions (|11|) of t, p, and 
q, but it is by no means the unique Hamiltonian that does this. With no 5^Pt specified a 
priori, the evolution of p t can take any form that we like, and this freedom translates into an 
ambiguity in both the Hamiltonian K\ and the energy E t . We explore the range of possible 
Hamiltonians (and therefore energies) by considering the functions that may be added to 
Kx without affecting the evolution equations (fTTjl. 

To start one would consider functions of all possible variables and parameters — that 
is functions of the form f(t,p t , q,p, X a ). However, we immediately see that a dependence on 
Pt, q, or p will change the evolution equations (fTT|) . Therefore only functions of the form 
f(t, A c ) may be considered. For such functions, the constraint equation K t = transforms 
to become 

p t + E t + ^f = 0. (18) 

Next, we demand that the constraint equations do not depend on the Lagrange multiplier 
A . This is equivalent to requiring that either the energy E t should not depend on A G or 
equally the equations of motion should not restrict the allowed values (or evolution) of A . 
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With this assumption, the freedom is reduced to f(X ,t) = \ g(t) + h(t), for any functions 
g(t) and h(t). Then, the derived equation of motion for p t ()15|) becomes 

^ = -4-1- < 19 > 

This will only be consistent with (fTHj) if h(t) is in fact a constant C . Then any Hamiltonian 
of the form 

*A = Ao(ft + |^ + ^ + <#)) + C , (20) 

will be consistent with our requirements, and so the energy will only be defined up to a free 
function: 

3=£ + X + »M- < 21 > 

We finish this analysis of the freedom using a physical argument. Mathematically, any 
function g(t) will satisfy our requirements. Physically however, it is reasonable to demand 
that this function should be in some way connected to the system. With that requirement 
we are reduced to considering functions of k(t). 



3. Gravity 

We now turn to a canonical Hamiltonian formulation of general relativity over a quasi- 
local region of space-time. Although the technical details will, of course, be much more 
complicated, many of the basic conceptual issues relating to extended phase space have 
already been dealt with in the previous example. Thus, for gravity we will also extend the 
usual phase space to include a time variable — which will be defined only on the boundary 
- and it will be joined by the spatial coordinates of the boundary. Further, just as k in 
the above was only fixed up to changes in the time parameter, for quasi-local gravity our 
boundary conditions will only be fixed up to intrinsic diffeomorphisms of that boundary. 
We will also find that a range of Hamiltonians will generate the standard evolutions of a 
space-time and that this freedom may be traced to the freedom to choose the evolution of the 
conjugate momenta of the boundary coordinates. Each of these Hamiltonians will be valid 
on its own constraint surface, and again each of these constraint surfaces will correspond to 
a different energy function for the system. 



3.1. Manifold without boundary 

Let us begin by briefly reviewing the Hamiltonian formulation of general relativity for a 
manifold with no boundaries. This will also serve to fix our notation and conventions. We 
will consider time-dependent fields living on a space-like three-manifold S. The intrinsic 
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geometry of E is fully specified by a space-like three-metric h a b] the derivative operator 
compatible with the metric will be denoted by D a . In the standard Hamiltonian formulation, 
the 3-metric h a b serves as the configuration variable, while its conjugate momentum is the 
tensor density P ab . Thus, the phase space consists of pairs of fields (h a b, P ab ) with the 
symplectic structure 

fi(<M 2 ) = / d 3 x {(5^(5^) - {5zh ab ){8 x P ah )} . (22) 

To specify a time evolution on E, we introduce a lapse function N and a shift vector 
field V a G TE. In the usual way the lapse and the shift will prescribe how time "flows" on 
E. Only after these fields are given can we define a time derivative 4 over E, as only then 
will we know how to associate points at a time t with points at time t + St, and also know 
how much proper time has passed during that interval. Given a flow of time (or equivalently 
a lapse and shift), we introduce a Hamiltonian 

H t = I d 3 x{NH + V a H a } (23) 
Jt, 

which generates time evolution. Specifically, we obtain 

^hab = [h a b](N,v) and (24) 

jP ab =[P a % Ny) , (25) 

where the exact forms of \h a b]{N,v) an d [P ab ](N,v) are given (along with expressions for the 
Hamiltonian constraint TC and diffeomorphism constraint TC a ) in |Appendix A Furthermore, 
the initial data for h a b and P ab must satisfy the constraints 

U = and H a = . (26) 

which are then automatically preserved in time. 



3.2. Boundary Conditions 

We would like to extend the Hamiltonian formulation of general relativity to manifolds with 
boundary. Therefore, we will now consider E to be a space-like 3-manifold with a closed 
2-boundary B. In this section, we will describe the boundary conditions enforced on B. The 
essential idea is to keep the boundary metric, lapse and shift fixed, up to diffeomorphisms 
of the boundary. To make this precise, we proceed as follows. 

Construct a three-manifold B = B x IR which is foliated by two-manifolds Bi = B where 
t is the foliation parameter and "=" indicates that the manifolds are diffeomorphic. The 
parameter i provides a notion of time on B. Specifically, the £>j are taken as "instants" of 
time, and i 2 is said to occur after t\, if i 2 > t\. We then introduce a time-like three-metric 7 a & 
on B which pulls back to a space-like two-metric a a b on each of the elements of the foliation. 
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Using the three-metric, we can also obtain the future-directed unit normal vector field to the 
Bt, which we denote u a . Finally, we introduce a time evolution vector field T a G TB which 
satisfies 

Cft=l. (27) 

We note that T a may be decomposed into its parts perpendicular and parallel to the Bi as 

f a = Nu a + V a , (28) 

where N is the lapse function and V a G TBi is the shift vector on B. Then, the future 
pointing unit normal to Bi is 

u a = -Ndi a . (29) 

We can also decompose the metric ^ a b as 

lab = -N 2 dt a dt b + o ah , (30) 

while 

-ab = -S-f a f b + -Lf (a V b) + a ab . (31) 
With this framework in place, we are ready to introduce our boundary condition on B. 
The Boundary Conditions 

(i) Construct the 3-manifold B and equip it with a foliation, fixed time-like boundary 
metric i a b and time evolution vector field T a as described above. 

(ii) Introduce a smooth diffeomorphism 

u : B x R -> B, (32) 

so that if t is the time parameter for S and B, and i(t) is some monotonically increasing 
function from IR — > H, then 

u(;t) : B ^ B m , (33) 

is a diffeomorphism for all t G IR (see figure |IJ). That is, instants of time on B are 
mapped to our pre-defined "instants" of time in B. 

(iii) The two-metric a a b, lapse N, and shift V a on B at a time t are equal to the corresponding 
a ab , N, and V a pulled-back to B using co(-,t) (or pushed-forward using cj _1 (-,t) in the 
case of V a ). 

(iv) The extrinsic curvature induced on B by /i a j and P ab must satisfy the time-like 
diffeomorphism constraint. 
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In condition (iii), we have required V a = V a which ensures V a G TB on B. However, in 
the initial set-up, we only required that V a G TS. This extra restriction is equivalent to the 
orthogonality assumption of Brown- York (which is discussed in more detail in section Q| and 
is not essential but will somewhat simplify the already involved discussion that will follow. 
The more general non-orthogonal case, which only requires that N 2 — V a V a = N 2 — V a V a 
and is essential if we wish to allow B to be null or space-like, has been studied and will 
be discussed in a future paper. For now, however, we will work with this orthogonality 
assumption and so find that 

l = £ " (34 > 

when these operators act on fields that are defined on B and then mapped back to B. Finally, 
(iv) says that even when we later consider general variations of our fields, they must continue 
to satisfy the diffeomorphism constraint at B. 

In order to describe the phase space of quasi-local general relativity, we need a concrete 
realization of the diffeomorphism uj. To obtain this, we provide a coordinatization of the 
three-manifold B in terms of the foliation parameter t and two "angular" coordinates x A on 
Bt §, which are chosen to be "time- independent" — that is 

C T x A = 0. (35) 

(Here, and in the remainder of the paper, we will use capital letters to signify the coordinates 
and components of tensors in this chosen coordinate system, while lower case letters are the 
usual abstract index notation.) Once the coordinate system on B is given, u may be described 

§ In many cases, multiple coordinate systems will be required to cover B without coordinate singularities 
and corresponding singularities in the coordinate vector fields. The calculations of this paper extend directly 
to cover these situations, but in the interests of clarity, we'll proceed as if one set of angular coordinates was 
sufficient. 
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by how it pulls back the coordinate system to B x JR. Giving (x A , t) at each point in B x 1R 
uniquely determines the map uj. Thus, we can (and do) specify uj by assigning an x A and t 
to each time and place in B. 

Notice that while (ii) specifies that a diffeomorphism exists, the exact form of uj is not 
fixed. However, any permissible uj can be obtained from any other || simply by composing 
uj with a diffeomorphism 

4>:B^B mapping (p, t) i-> {p', F) = <f>(p, t) , (36) 

which preserves the foliation of B. Thus, the freedom in defining uj is the freedom to consider 
a preferred uj composed with all 0, where 

o uj : B x R, — > B maps (p, i) i— > (a> (p, i)) . (37) 

Since u; D mapped all fields (t7 a f,, JV, K a , and the coordinates t and x A ) back to 23, then <pouj 
will too. That is, from a computational point of view, our diffeomorphisms act so that 
observers fixed to points either on B or B will see all of these fields perturbed. 

For an infinitesimal diffeomorphism 0, it is straightforward to calculate the changes to 
the various fields at the boundary. We begin by recalling that infinitesimal diffeomorphisms 
are generated by non-singular vector fields and so to each we may associate a vector field 
8X = (5i,5x a ). The requirement that preserve the foliation means that 8t should be 
a constant on each leaf of the foliation, although it can vary from one leaf to the next. 
Meanwhile, Sx a can be any non-singular vector field in TB that is everywhere parallel to 
the foliation surfaces. Under the infinitesimal action of this diffeomorphism the fields on the 
boundary B will change according to 

lab >-»• lab + £ SXlab 

dt I— > dt + C sxdi 
f ^T + C 5X T 
t i— > i + 5i and 
x A x A + 5x A . 

Here, 5x a is simply the vector field 5x, expressed using the abstract index notation, while the 
5x A are the components of the vector field 5x. It will be important to keep this distinction 
clear. The two are related by 

Sx A = (5x) a d a x A , (38) 

where d a is the intrinsic derivative over B. 

| Up to potential global topological obstructions which we shall ignore with impunity since in our calculations 
we will only consider "infinitesimal" variations of uj. 
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The corresponding change to the fields in B x H is given by 

N ^ N + (5t)^N + C 5x N (39) 

V* ^V a + (5t)^V a + £ sxV a (40) 
at 

o~ab o~ab + (St)—a ab + C Sx a ab , (41) 

t h^t + 5t and (42) 

x A i-> x A + 5x A , (43) 

after uj maps everything back to B. Note that we have dropped the over-bars in order to 
reduce notational clutter. In the future we will blur the distinction between SxR and B. 
Whenever there is an ambiguity, u is understood to be acting to make the identification and 
map quantities back and forth. 

To summarize, our boundary conditions imply that once a map u from B x 1R to B 
is specified, the boundary metric, lapse and shift are known. We have simplified matters 
by introducing a coordinate system on B which allows us to easily characterize the map uj, 
but the results obtained in the following subsections will not be sensitive to this coordinate 
system, and it is likely that the calculations could be done without introducing coordinates 
at all. What will be important is that the only allowed variations of the boundary will be 
infinitesimal diffeomorphisms. They will generate changes in the boundary metric, lapse 
and shift as given by (j39H41j) . Additionally, they will change the map u, or equivalently the 
coordinates associated to the points of B, according to 



Comparing with earlier work, Regge and Teitelboim |2j showed that their Hamiltonian 
was covariant with respect to the Poincare group at infinity acting on the boundary. Here we 
have set things up so that we may study the effect of the diffeomorphism group that maps 
the boundary into itself (while preserving the foliation) on the Hamiltonian. Thus, apart 
from the difference between boundaries at finite difference and infinity, we are also studying 
different group actions on those boundaries. 



3.3. Phase Space and Hamiltonian Evolution 

In the previous subsections, we considered the fields in the bulk as well as the boundary 
conditions imposed at B. In this subsection, we will turn our attention to the phase space 
and associated symplectic structure, as well as the Hamiltonians associated with evolution 
equations. 

With the boundary conditions that we have imposed, the configuration variables of a 
point in the phase space will be given by the three-metric h ab , a value of the time parameter 
t, and a set of angular coordinates x A on B. The conjugate momenta to these configuration 
variables will be P ab , P t , and Pa, and so the coordinates of a point in phase space will be 
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given by the six fields (h a b, P ab ,t, P t , x A , Pa)- Note that our boundary conditions say that 
once we specify t and the x A , we will know a a b over B, thus h a b cannot be chosen completely 
independently of those variables. The symplectic structure % on the phase space then takes 
the form: 

^(Si,6 2 ) = J s d 3 x{(6 1 h ab )(5 2 P ab ) - (6 2 h ab )(6 1 P ab )} (44) 

+ J s d 2 x{(6 1 x A )(6 2 P A ) - (6^(6^)} 
+ {(Sit)(S 2 P t ) - (6 2 t){6 1 P t )} . 

The variations which appear in the bulk are entirely free, and are not restricted to being 
on-shell. However, they are partially restricted at the boundary Specifically, variations must 
preserve the boundary conditions, and so must be generated by some (St, Sx A ) . Then the 
variations of the lapse N , shift V a , two-metric a a b, the time parameter t, and the coordinates 
x A are given by equations f)39H43j) . 

A typical on-shell variation S\ will be generated by the action of 



A = A - + £ A , (45) 
at 

where A Q must be constant on the boundary (so as to preserve the foliation of the boundary) 
and A G TH t in the bulk while it is tangent to B t on the boundary. Then we know how our 
phase space variables must evolve with A. Namely, 

SAhafi := Xo — hab + £ x hab, (46) 

5 k P ab ■= \ j t P ab + £ \P a \ (47) 
S A t := A , and (48) 
S A x A :=\ a d a x A , (49) 



where the time derivative of h a b and P ab is given in ()24j) and (J25j) respectively. The evolution 
of the remaining two fields, Pt and Pa is currently undetermined. 

We would now like to determine whether there is a Hamiltonian K\ which generates 
the on-shell evolution S A given in equations (|45H49|) . Thus we will try to manipulate the 
symplectic structure into the form 

fl(S\, S) = 8K\ + constraints. 

In doing this we will be aided by the fact that 5^P t and S^Pa are as yet unspecified. 

This is really the pre-symplectic structure since the existence of the Hamiltonian and diffeomorphism 
constraints means that we have not properly isolated the true degrees of freedom. The proper symplectic 
structure would restrict itself to these real degrees of freedom. That said, we will follow the standard practice 
and work with the pre-symplectic structure. 
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The first important step in this process is to evaluate the bulk term in the symplectic 
structure. Here, we make use of the result of | Appendix B 



, d 3 x { (5 A h ab ) (6P ab ) - (5 A P ab ) (6h ab ) } = (50) 
5 (J d 3 x {\ NH+ (X a + \ V a )H a } + J d 2 x^ [X (Ne - V a Ja ) - X a Ja 

- / d 3 x{6(\ N)H + 6(\ a + \ V a )H a } 

JT, 

— \ d 2 x\fa 

JB 



eS(\ N) - Ja S(\ V a + A") 



XnN 



s da, 



As one would expect from its general appearance, the derivation of this result is fairly 
involved. It is discussed in greater detail in the appendix. Here however, we simply note 
that 



3a 
s ab 

-ab 



k/(8irG), 

- 2a ac P cd n d /\/h, and 
(1/8ttG) (k ab - (k - n c a c )a oh ) 



(51) 



where k = —o~ D a n b is the extrinsic curvature of B E £ (note the sign convention and recall 
that D a is the covariant derivative in S that is compatible with h ab ), a c = jjD c N and G is 
Newton's constant. The meanings of some of these quantities are easier to see if one thinks 
of them as being defined in four- dimensional space-time — this perspective is discussed in 
section H| following equation (JH3J. 

In the last line of (|5U|) . there are terms involving the variation of the boundary 2-metric, 
lapse and shift. However, due to our boundary conditions, we know that 5 must act on the 
lapse, shift, and two-metric like an infinitesimal diffeomorphism in B generated by the vector 
field 

5X a = 5tT a + (5x) a . 

Then, with the help of |Appendix C[ we see that because the equations of motion hold at the 
boundary, 



(52) 



f B d 2 x (sC sxN - j a C sx V a - — s ab C sx a ab ^ 

I d 2 x(6t)^-(^[Ne-V a j a ])- I <Px(6x)*±{y/tj a ) 
JB t dt JB t dt 



where the four-dimensional, in-boundary Lie derivatives C t from |Appendix C are replaced 
with from our three-dimensional perspective. We can now use ()50|) and (|52[) to rewrite 
the bulk part of the symplectic structure in (J44)) to obtain 

Q(8 A , 5) = 5 (J d 3 x [(X N)H + (X V a + \ a )H a ) 
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+ 6 U d 2 x [Xo(Ne - V a ]a ) - X a j a ] + X Q P t + J al 2 x X A P A 

- [ d 3 x [5(X N)H + S(X a V a + X a )H a } 

- (5X ) \p t + / d 2 x^(Ne - V a 3a ) 

JB 

-{St) S A P t + I d 2 xX ^-(^(Ne-V a 3a ) 
jb at v 

(5X A )P A -(5X a )V^Ja 



d x 



d 2 x (5x) 



5 A P A -X — (V^Ja) 



(53) 



This is nearly the desired form. We see that the first two lines are an exact variation, as 
desired. The third line will vanish provided the usual bulk constraints of general relativity 
are satisfied. The fourth line is then a new constraint which appears at the boundary and 
relates P t to other fields on the boundary. The fifth line gives the evolution of P t and 
guarantees that the constraint will continue to hold (these two correspond to the constraint 
and evolution equation found for p t in the case of the harmonic oscillator) . We would like to 
write the final two lines in a similar form — a constraint for P A and an evolution equation 
which ensures that this constraint is preserved. To this end, we must evaluate S(X A ). Making 
use of (13~%1). it follows that 



5(\ A ) = (5X) a {d a x A ) - X a d a {5x A ) 



(54) 



where as usual the capital Latin index indicates a component of a field in the coordinate 
system imposed on B, while the lower case index is an abstract tensor index. Therefore, we 
obtain: 



S(X A )P A - (SX) a V^ja = (Pa - V^3a)S(X a ) - ^ ]A X a d a (S 



x 



(55) 



where j A = j a ^ 
parts, and using Sto 



8 



Substituting this expression into (jB^j) . doing some integration by 
£es theorem, we can rewrite the last two lines as 



- / d 2 x(5X A ) \P A - V^Ja] - / d 2 x (Sx) A \S A P A - 5 A (^ JA 



(56) 



where as usual 5 a = A ^ + £ \. 



In order to simplify matters, let us introduce some notation: 

h a = f d 3 x {x (nh + v a n a ) + x a n a } 

Jt, 

K t =P t + [ d 2 x^/a[Ne - V a j a ] and 

JB 

La — Pa ~ \f^]A ■ 



(57) 

(58) 
(59) 
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Making use of the new notation, as well as (|56jh we can rewrite the symplectic structure 
(JS3J) as 

Q(5 A , 5) =6 (h a + \ K t + J d 2 x \ A L A ^j 

- [ d 3 x [5(\ N)H + 5(X V a + \ a )H a ] 

- (SX )K t - [ d 2 x(5X A )L A 

JB 

- (5t)(5 A K t ) - / d 2 x(6x) A (6 A L A ) . (60) 

JB 

In order for a Hamiltonian generating evolution along A to exist, the right hand side of (jHOj) 
must be an exact variation. Thus, all terms after the first — which is already an exact 
variation — must vanish. This can be accomplished if the last three lines vanish. Therefore, 
the evolution along A will be generated by a Hamiltonian provided that we are "on-shell," 
by which we mean: 

H = and H a = 0; 

K t = and L A = . (61) 

The first two expressions are the usual constraints of general relativity, while the last 
two restrict our evolution to a constraint surface in the extended phase space. They can 
equivalently be thought of as fixing the values of P t and P A to be 

Pt = ~ I d 2 x ^[Ne - V a Ja ] and P A = V^Ja ■ (62) 

JB 

Finally, the last two terms in (|5U|) define the action of 5 A on P t and P A . Essentially, it fixes 
the evolution so that it evolves points on the constraint surface defined by equations (|H2*|) 
into other points on that surface. 

Therefore, we have shown that "on-shell," i.e. when (J6T| is satisfied, the evolution along 
the vector field A is generated by the Hamiltonian K A which is given as: 

K A = [ d 3 x{(\ N)H + (\ a + \ V a )H a } + \ K t + [ d 2 x\ A L A . (63) 

JS JB 

Although we have introduced a coordinate system on the boundary in order to characterize 
the diffeomorphism lu, the final form of the Hamiltonian is independent of this choice. 
Specifically, the term X A L A will be the same when evaluated in any set of coordinates. 
As in the previous example, this Hamiltonian is not determined uniquely. However, we shall 
postpone discussion of the ambiguities to the next subsection. 

For the spherically symmetric similar analysis to this one may be found in [Jj. 

In that case, the symmetry means that there is no need to include "angular" coordinates as 
configuration variables. 
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3-4- Energy and Angular Momentum 

The on-shell value of the Hamiltonian ()63)) generating evolution along A will be zero. To 
find the energy and angular momentum we must again consider the value of the conjugate 
momenta rather than the Hamiltonian. Since the equations of motion guarantee that K t 
vanishes on-shell, we can use (I58J) to find the value of Pt. The energy associated to time 
translation is then simply the negative of this: 

E d/dt :=-P t = J d 2 x (Ne + V a Ja ) , (64) 
Similarly, the equations of motion show that La vanishes on shell, so that 

Pa = V^Ja, (65) 
Combining these two results, we arrive at the expression 

E A = f d 2 x {X (Ne + V a Ja ) + X a Ja } , (66) 

JB 

for the "charge" associated with evolution according to A. Specializing to the case of pure 
time evolution, we obtain (jHlj) which is the energy associated with the time translation 4. 
Similarly, given a vector field <j) tangent to the 3-surface S and tangent to B at the boundary, 
the angular momentum is given by 

J = / d 2 x^ a Ja ■ (67) 

JB 

As in the case of the Hamiltonian, our expressions for the energy and angular momentum 
associated to the boundary were obtained by making use of a specific coordinatization of the 
boundary. However, the final results are independent of the choice of coordinates, as one 
might have hoped. 

Finally, we turn to the ambiguity in the Hamiltonian K\ and the corresponding freedom 
in the definition of the energy and angular momentum associated to the boundary. We 
follow the same procedure that we used in section 12.21 to explore the allowed freedom in the 
Hamiltonian for the simple harmonic oscillator. Namely, we start by considering adding any 
functional of the form 

F = [ d 2 xf(h ab , P a \ t, P t , x A , P A , A , A a ) , (68) 

J B 

to K\. However, on requiring that: 

(i) the new Hamiltonian still generates the same evolution equations 1)461149)1 . 

(ii) the new versions of the constraints Kt — and La = remain independent of the 
Lagrange multipliers A Q and X A (equivalently either energy and angular momentum are 
independent of these parameters or A G and X A may be freely chosen without reference 
to the equations of motion or constraints), and 
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(iii) the new versions of these constraints are preserved by the new evolution equations for 
P t and P A , 

the freedom is greatly reduced. It turns out that subject to these conditions, the only 
valid Hamiltonian functionals will take the form, 

K' A = K A + X F(t) + C , (69) 

where F(t) is a free function and C is a free constant. 

Furthermore, if we again make the argument that the only physically relevant functions 
F(t) are those which are associated with the boundary data, then these free functions must 
take the form 

F(t)= f d 2 x^f(a ab ,N,V a ), (70) 

JB 

where / is a free function of the given fields. Therefore, we find that the energy associated 
to an evolution along A is 

E' A = / d 2 x [X (Ne + V a j a ) + X a j a ] + / d 2 x A oV ^/(a afe , N, V a ) . (71) 

JB JB 

Readers familiar with will immediately recognize Ed/dt with this F(t) as the Brown- York 
energy, including the usual reference term. In contrast, there is no freedom to add a reference 
term to the angular momentum expression. Its most general form is 

4 = J B d 2 x^ a Ja, (72) 

unchanged from equation (|67j). In the Brown- York derivation a freedom remains in this 
definition. 

Although specific choices of F(t) will not be our major concern here, a few words are 
in order. The F(t) is a term in the expression for energy that cannot be determined by the 
formalism. Indeed as far as our formalism is concerned all F(t) are equally valid. Thus, 
a particular F(t) can be chosen to suit the needs of a particular problem. Popular choices 
in the literature are usually based on the requirement that Ed/dt vanish for some particular 
space-time (for example flat space or AdS). The reader is directed to any of the papers listed 
in [Sj for further discussion of this point and other applications of this energy expression. 

4. Comparison with action arguments 

Let us now compare our calculations with those used in the Brown- York derivation of the 
quasi- local energy. Their original paper [S| pursued two lines of argument, both of which 
follow from an analysis of the gravitational action / and the boundary conditions which must 
be imposed on quasi-local boundaries so that the action principal will be well defined. From 
the action they proceeded in two directions to obtain the quasi-local energy. The first used 
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Hamilton- Jacobi type arguments and followed a more qualitative path, while the second used 
a Legendre transform to derive a Hamiltonian functional from the action. In the next few 
paragraphs we will review the Legendre transform route. 

To begin, let us recast the three-dimensional constructions of the previous section into 
four- dimensional space-time. We consider a region of space-time M that is bounded by a 
time-like boundary B and two space-like boundaries £j and S 2 as shown in figure EJ The 
metric over M will be g a b and the associated covariant derivative is V a . The £ of the last 
section is now replaced with space-like surfaces E f that foliate M. We assume that Si and 
£2 are leaves of that foliation. The S t induce a foliation Bt of B, and these Bt replace B. 
The boundary metrics h a b, jab, and a a b are all induced by g a b and may be calculated as 

Kb = 9ab + u a u b , (73) 
lab = g a b - n a n b , and (74) 

(Tab = gab + U a U b - U a n b , (75) 

where n a is the outward-pointing space-like unit normal to B and u a is the future pointing 
time-like unit normal to the Sf. In defining these metrics we have built in the orthogonality 
assumption that u a n a = 0. Allowing non-orthogonal intersections increases the complexity 
of the calculations but doesn't substantially change the following results [12 . 

To understand the Brown- York arguments, let us begin with an analysis of the well- 
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known trace-K action for the quasi-local region M: 

I = — [ d^x^g'K + - / d 3 xVhK - - [ dPx^jQ . (76) 
2k Jm k Jt, k Jb 

In the above, TZ is the Ricci scalar corresponding to g ab . Following the Brown- York sign 

convention for extrinsic curvatures, K = —V a u a is the trace of the extrinsic curvature of 

Si.2 in M while O = — V a n a is the trace of the extrinsic curvature of B in M. Further 

is = Je 2 - Jsx and K = 87rG - 

The first variation of this action with respect to the metric g ab is (see any standard text 

on relativity, for example [3]), 

SI = ^- [ d 4 x^g-(G ab +Ag ab )6g ab + f d 3 x (P ab 6h ab ) + f d 3 x (it ab 5 lab ) .(77) 

Note that in four-dimensions, P ab is not an independent field as in the Hamiltonian 
formulation, but is instead 

P ab ■= 2p [Rh ab - K ab ) , (78) 

where K ab = —h a c h b d V c Ud is the extrinsic curvature of Si j2 and K is the trace of K ab as 
discussed above. Similarly, 

7r ab ■= -^El (e 7 afe - e a6 ) (79) 

is an equivalent tensor density defined by the surface B, with Q ab = — 7 a c 7b d V c nd being the 
extrinsic curvature of B. 

Then, if the boundary metrics h ab and ^ ab are fixed so that 

5j ab = and 5h ab = 0, (80) 

the variation of the action vanishes on-shell. Note however, that this means that rotations 
or translations of the boundary are not allowed variations of the action, unless those 
rotations/translations are generated by Killing vector fields. We will come back to this point 
momentarily, but for now we note that requiring that the variation of the action vanish and 
boundary conditions be met gives the Einstein equations in the standard way. Further note 
that with these boundary conditions, free functionals of the boundary metrics may be added 
to the action without affecting its first variation — since those terms are fixed, the variations 
of the functionals will vanish. As such, the exact form of the action is ambiguous and 

r = i + i { lab ,h ab ], (si) 

is an equally valid action for any functional I Q of the boundary metrics, in the sense that 
fixing the same boundary conditions and demanding that 51' = will also generate the usual 
equations of motion. 
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Next, in preparation for applying a Legendre transform and so obtaining a Hamiltonian, 
we take note of the foliation of M. This gives an (imposed) notion of "instants of 
simultaneity." A time flow is put on M in the guise of a vector field T a that satisfies 
T a d a t = 1 and which lies in TB everywhere on B — that is, the time flow generates the 
boundary. From our orthogonality assumption that u a n a = 0, we know that u a is both the 
unit normal to the S 4 as well as to the B t in B. Thus, we may write 

T a = Nu a + V a , (82) 

for some lapse N and shift vector field V a that is everywhere an element of TBt over B. 
Making use of this time evolution, we can perform a Legendre transform on the action to 
obtain: 



where 



/ = Jdt^-Ht + J^ al 3 x (P ah C T h ab ) | (83) 
H t = [ d 3 x\NH + V a H a }+ [ d 2 x^/a(Ne -V a j a ). (84) 

JSt JB t 

Now, e retains its earlier value (see equation ([51))) but in the context of this four- dimensional 
approach j a can be written as 

3* = ^-^ a o. b u c V b n c , (85) 

07TCt 

and so is proportional to the connection on the normal bundle to B t . Referring back to (|51|) 
it is probably also worthwhile to note that in four-dimensions a a = u b V b u a is the acceleration 
of the unit normal vector field u a along its length. 

The quasi-local gravitational Hamiltonian for the region of space-time M is given by 
H t . However, there is no reason that the functional I [lab,hab] appearing in (J8T|) should 
decompose into an integral with respect to t and so give us a Hamiltonian reference term 
that is local in time — extra assumptions have to be made in order for this to be true. The 
usual requirement is that I Q be a linear functional of the lapse and shift [U] which guarantees 
that the energy density e and angular momentum density j a are independent of lapse and 
shift. Then functions of the form 

H = / d 2 x {Nf(a ab ) + V a f a (a ab )} , (86) 

JB 

where / is a free function of a ab and /„ is a free vector valued function of the same, may be 
added to the Hamiltonian. 

Comparison of and shows that the phase space methods of section El and the 
action arguments given here derive the same quasi-local energy, up to assumptions made 
about the boundary conditions. Recall that in the phase space calculation, we found free 
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functionals of time only — demanding a connection with the boundary data allowed us to 
write them in the form of equation (|70j) 

F(t)= f d 2 xf[a abl N,V a ]. (87) 



By contrast, the Brown- York freedom is more general and allows free functionals of 
the boundary metrics in the action without any need to integrate. It is only after 
extra assumptions that they can be broken up into Hamiltonian reference terms that are 
independent of local details of the geometry of B. This difference in boundary terms arises 
because the boundary conditions that we have imposed are not the same. The Brown- York 
conditions are more restrictive than ours (and so allow more freedom in the free functionals). 

To obtain a fairer comparison, let us weaken the boundary conditions in the action 
formulation so that they are equivalent to those used in section EJ 

(i) Instead of rigidly fixing the metrics on the boundaries, we allow variations which are 
foliation preserving translations/rotations of the boundary. Thus, the overall geometry 
of the boundaries will be fixed, although particular features of the geometry will not be 
fixed to particular points of the manifold. To do this, we allow the variations to act as 
diffeomorphisms on the boundary metrics, with the restriction that they should map B 
into itself. Then, 



S'jab = C Ylab and 6h ab = £ z h ab , 

for some vectors Y a £ TB and Z a £ TS for which Y a = Z a £ on those corner 
two-surfaces. 

(ii) At the same time, we impose the new condition that the diffeomorphism constraint 
should hold on the boundary surfaces Si, S 2 , and B. That is 

D a P ab = and A a 7r afe = 0, (89) 

where, as before, D a and A a are the induced covariant derivatives on Ei^ and B 
respectively. Of course, these constraints hold automatically on-shell. In the action 
formulation however, variations are not restricted to being on-shell and we will need to 
impose these conditions in the following calculation. 

Then, with these new conditions: 

61= f d 3 x7i ab £ Ylab + f d 3 xP ab £ z h ab (90) 

JB JSo-Si 



2 / d 3 xir ab A a Y b + [ d 3 xP ab D a Z b 

JB 

2 f d 3 xA a (n ab Y b ) + f d 3 xD a (P ab Z b 

JB 



d 2 X\l <~r 

B2-B1 



,u a 7T ab Y b + l=n a P ab Z b 
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1 



( 



L 



al 2 xy/aY b [u a W b n a + n a W h u a ] 



B 2 ~B 1 



0. 



The second line applies the representation of a Lie derivative in terms of the covariant 
derivatives A a and D a of B and £^2 respectively. The third uses the fact that the 
diffeomorphism constraint holds on the boundary. The fourth integrates bulk terms out 
to the boundaries of B and £1,2, while the fifth uses the fact that Y a = Z a G TB12 on 
those boundaries. Thus, with these modified boundary conditions, the action also vanishes 
on-shell. + 

Given these boundary conditions, there are other actions for which 51 = and also give 
the same equations of motion. In particular, consider any reference term which is of the 
form 



where t is some labeling of the foliation which is equal to ti on £1 and £2 on £2. F(t) is a 
functional that depends on the foliation surface only. Then SI = for the variations that 
we have considered. On applying the Legendre transform, we will find that for I + I Q the 
Hamiltonian becomes 



which is the same form that we found with our phase space analysis. Thus, when we 
impose the same boundary conditions for the two approaches, we obtain the same value 
for the Hamiltonian functional — as would be expected. As we have seen above, with the 
standard Brown- York conditions which fix the metric at each point in the boundary the 
free functionals can depend on specific local features of the intrinsic geometry (equivalently, 
we could think of a coordinate system imposed on the boundary with the functional being 
allowed to depend on the coordinate system as well as the metric). With the new boundary 
conditions however, that freedom has changed and the free functionals can only depend on 
global geometric features (for example integrals of geometric invariants over the boundary). 

5. Discussion 

In this paper we have introduced a phase space description of general relativity in a manifold 
with boundary. The boundary is equipped with a preferred foliation and a metric which 
is fixed up to diffeomorphisms which preserve this foliation. Therefore, we must allow 

+ The essentials of this proof that the action is invariant with respect to diffeomorphisms that map the 
boundary into itself may be found in |14|. where they occur in the course of a demonstration that the 
Bianchi (and other) identities may be derived from the action principal. 




(91) 



H t = F(t) + [ al 3 x [NH + V a H a ] + [ al 2 x^{Ne - V a j a ) , 



(92) 
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phase space variations which act as diffeomorphisms of the boundary metric. This can 
be accomplished by introducing a coordinatization of the boundary and extending the phase 
space to include both the time and spatial coordinates of the boundary and their conjugate 
momenta. Then, in the extended phase space, we are able to discuss variations which induce 
diffeomorphisms of the boundary, even if they are not symmetries of the boundary metric. 
In particular, time evolution is an allowed variation even if the boundary metric is time 
dependent. Similarly, rotation along a vector field which is not a Killing vector of the 
boundary metric is also permitted. 

Furthermore, we are able to associate "conserved charges" with these motions. In the 
case where the motion is not a symmetry of the boundary, the "conserved charge" will not 
be constant. Hence, the conserved charges cannot be equal to the on-shell value of the 
Hamiltonian (which is by definition a constant). Instead, the on-shell values of the momenta 
conjugate to the coordinates give the "conserved charges" associated with these motions. In 
particular, if the variation induces a time translation at the boundary, the corresponding 
conjugate momentum is equal to the energy of the boundary. Similarly if the variation 
produces a rotation of the boundary, the conjugate momentum is the angular momentum 
of the boundary. These quantities need not be conserved: they can vary as you move up 
the boundary. In addition, the final expressions for the energy and angular momentum are 
independent of the choice of coordinate system of the boundary, as one would hope. 

The expressions for energy and angular momentum obtained from the Hamiltonian 
formulation are the same as those obtained in the Brown- York formalism (up to the free 
functionals where a difference arises due to the slightly different boundary conditions). This 
is unsurprising since both methods describe the same situation, just in a different language. 
In particular, the energy obtained is not uniquely fixed and there is the freedom to add to 
it a function of the boundary data. When we harmonize boundary conditions as we did in 
the last section, this is the same freedom as found in the Brown- York formalism and can be 
fixed in the usual manner. It should be possible to extend the formalism presented here to 
the case of non-orthogonal boundaries as considered in > and work is already under way 
to do so. * Additionally, we have required throughout that variations preserve the foliation 
of space-time. It is likely that this condition can also be weakened or removed entirely. 

The approach we have described above should have several applications. Indeed, it was 
initially motivated by a desire to extend the isolated horizon framework to include the 
physically interesting case of dynamical horizons. In order to describe these horizons within 
a Hamiltonian framework, it is necessary to allow for a varying black hole mass and angular 
momentum. Making use of the extended phase space methods introduced here will make 
this possible ^B] . A second application would be to consider the Hamiltonian formulation of 

* In fact, in his alternative formulation, Kijowski [3] has already shown that the non-orthogonal "angle" 
parameter may be included as a configuration variable in the phase space. 
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general relativity with a boundary at null infinity. Again, the mass and angular momentum 
are not constant as they can be radiated away. Therefore, the framework here would again 
be relevant. It would be interesting to recast the results of Wald and Zoupas [T2j in terms 
of the phase space description presented here. 
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Appendix A. Constraints and time evolutions 



For completeness we include the expressions for [h a b\(N,v), [P ab ](N,v), T~t, and TLb (in the 
absence of matter and a cosmological constant). Keep in mind that the tensor (and tensor 
density) fields in these expressions all live on £ and its assorted tensor bundles. Further, 
D a is the intrinsic covariant derivative over £ and k = 8tiG, where G is the gravitational 
constant. 



[h 



AkN 



ab\(N,V) 
pab 

i (N,V) 



ab 



;Ph ab ]+C V h 



ab ! 



(N ^G ab 
2k V 



D a D b N - h ab D c D c N )+C v P 



, ab 



•tab 



+ ^1 f[P cd P cd - -P 2 ]h ab - A[P< a P c V - -PP ab 



H 



Vh V 

2k V 
2D b P a b . 



P ab P ab - -P'< 



and 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



In the above ^G a b = Rab — \Rg a b-, where R a b is the Ricci tensor on S and R is its contraction 
with the metric h a b. 

Full derivations of the above results from the four-dimensional geometry of solutions to 
the Einstein equations can be found in most standard relativity texts, or using this notation 



Canonical Phase Space Formulation of Quasi-local General Relativity 



27 



in Essentially though, they come from thinking of E as one leaf in the foliation of 

a solution of the Einstein equations (M, g a b) into space-like hypersurfaces (as discussed in 
section EJ) . If u a is the future pointing unit normal to E, the time evolution given by (N, V a ) 
corresponds to that generated by the vector field T a = Nu a + V a . Further P ab is no longer 
an independent variable, but instead is closely connected to the extrinsic curvature of E in 
M (see equation (|7Sjl). 

The Einstein equations tell us that in empty space, G a b = 0. Then, the Hamiltonian 
constraint ()A.3jl is equivalent to the statement G a bU a u b = 0, the diffeomorphism constraint 
(jA.4j) is equivalent to h b a Gb c u c = 0, and the evolution equation for P ab (jA.2|) is equal to 
C T P ab with h c a h d b G cd = being used to rewrite the result in terms of quantities defined 
entirely on E. Finally, the evolution equation for h a b is equivalent to Crhab- 



Appendix B. Variational calculation 



In this appendix we show that if A D is any function on E which is constant on B and A a G TE 
is any vector field that lies in TB on B so that 

6 A P ab = X [P ab } {Ny) + C x P ab and 8 K h ah = \ [h ab ] (N , V ) + C x h ab , (B.l) 



then 



d 3 x 



(5 A h a b)(5P ab )-(5 A P a b)(5h ab ) 



(B.2) 



S (J afx [X NH + (A a + \ V a )H a } + J^al 2 x^ [\ {Ne - V a Ja ) - \ a Ja 
- J d 3 x [5(\ N)H + 5(X a + \ V a )H a ] 



d 2 x yfa 



e6(\ N)- Ja 6(\ V a + \ a ) 



—p^s da ab 



where 



£ 

ja 
„ab 



fc/(87rG), 

- 2a ac P cd n d /Vh, and 
(1/8ttG) (k ab - {k - n c a c )a a 



k a b = —° r a c o~b d D c n d is the extrinsic curvature of B in E (note the sign convention), k = a ab k ab , 



and a c = j^D c N . 



The calculations needed to obtain this result are quite lengthy, but luckily the bulk of 
them may be found in previous papers. The key result may be found in JH] or ^3] , and says 
that for a general variation S (there is no restriction to it being on-shell), 



5 



J^d 3 x[NH + V a H a 



(B.3) 
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J d 3 x {H5N + H a 5V a + [h ab ) {N , v) 6P ab - [P a \ Ny) 5h ab \ 
- J^d 2 x \n5{^s) - V a 8{^j a ) + \^Ns ah 8a a 



>ab 

Now, a careful examination of the calculations found in the above sources, shows that if we 
replace the N and V a on the left hand side of the equation with more general terms — for 
example any scalar field a over £ and any vector field (3 a over S that is restricted to be 
parallel to B on that boundary, then the expression changes to become 
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J^d 3 x[aH + (3 a H a ^j = J^d 3 x{H8a + H a 8f3 a } (B.4) 

+ / d 3 x{^h ab )(8P ab )~^P ab )(5h ab )} 
- J^d 2 x \a5{^e) - (3 a 5(y/aj a ) + ^y^as ab 5a ab j , 

where f = (a/N)f t + £ 

p-{a/N)v an d 8^P ab and 8^h ab are defined in an analogous way to 
5\P ab and S\h ab in equation (|B.1[) . Therefore, in order to obtain the desired result (jB.2|) . we 
must take 

a = X Q N , (3 a = \ a + \ V a , 
and re-express the surface term as 

J^d 2 x {a5(v^r) - (3 a 8{^j a ) + -^as ab 8(T ah \ = 

5 (J d 2 xy^ [as - p a j a )}\ - J d 2 x^ \e8{a) - 3a 5([3 a ) - ^s ab 8a ab 
to obtain (IB.2I). the desired result. 



Appendix C. The diffeomorphism constraint on the boundary 

In this appendix we consider a solution to the full four-dimensional Einstein equations over 
a region M like that discussed at the beginning of section We show that if 

(i) the diffeomorphism constraint holds on a time-like surface B that has intrinsic metric 

"fahi 

(ii) B is foliatied with closed, space-like, two-surfaces B t and there exists a vector field T a 
such that T a [dt] a = 1, and 

(iii) X a is a vector field taking the form X a = X T a + X a where X Q is a function of t alone, 
while X a is any non-singular vector field over B that is everywhere parallel to the Bt 
surfaces, 



Canonical Phase Space Formulation of Quasi-local General Relativity 



29 



then 

J b d 2 x [eC X N - Ja C x V a - ^s ab C x a ab 

J b d 2 xX C T (V^[Ne - V a Ja ]) - j B d 2 xX a C T (^j a ) . (C.l) 

where N and V a are the usual lapse and shift defined so that if u a is the forward-pointing 
unit normal to the B t (that is T a u a < 0), then T a = Nu a + V a . a ab is the induced two-metric 
on the B t surfaces, and e, j a and s ab are defined in (|51|). 

We begin by introducing the conjugate variable to ^ abl namely it ab which is given by: 

= M ~ ^ ' (C - 2) 
where Q ab = 7 a c 7fe d V c nd is the extrinsic curvature of B with respect to the unit normal n a . 
Then, we can write e, j a and s ab in terms of 7T ab SO that 

£ = - 2vr a Wv /3 7 7 
j a = 2a ab ir bc u c / , and 
s ab = - 2a a cr h 7r /V-7 • 
Now, taking A a as the covariant derivative on B we first note that the diffeomorphism 
constraint implies that 

A a n ab = 0, 

and therefore 

/ d 2 xA a (-K ab X b ) = ~ / rf 2 xvr ab /: X 7 ab . (C.3) 

The result (jQH) arises when we break up each side of the expression ()C.3|) into parts parallel 
and perpendicular to the foliation. To that end note that 

lab = o- ab - u a u h and v / = 7 = N^fa . 

Then, working on the left hand side we can show that 

V^A a (n ab X b /^) = -C T (V^u a n ab X b /^) 

+ £ v (V^u a 7r ab X b /^) + V^d a (Na a b 7r bc X c /^) . (C.4) 

Since the surfaces Bt are closed, when we integrate this expression only the first term on the 
right hand side survives, giving 

/ d 2 x A a (7i ab X b ) = \ f d 2 xC T (v^(X D [Ne - V a 3a ] - X a ]a )) . (C.5) 

J Bt 2 JBt v ' 

Let us now turn our attention to the right-hand side of equation (jC.3|) . We begin by 
recalling that u a = —N[dt} a . Therefore, with X Q constant on each slice and C x(dt a ) = 
(CTX°)dt a it is simple to show that 

C-xUa = ( ~ +C T X^j u a . (C.6) 
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In particular, it then follows immediately that o C x u a = 0. Finally, by decomposing the 
metric ^ ab in terms of u a and a ab we obtain: 

7T ab C xlab = v/=7 [e u a C x u a + Ja (C u X a + a ab C x u b ) + ^s ab C x a ah } 

= sfd [eC X N - j a C x V a - ^s ab C x a ab ^ - y^j a £ T X a 

+ ^(Ne-V a j a )C T X . (C.7) 
Equations ()C.3|) . ()C.5|) . and ()C.7|) can then be combined to give the promised result: 

jf al 2 x [eC X N - j a C x V a - ^s ab C x a ab ^j = 

jf d 2 xX C T (^\Ns - V a j a }) - jf d 2 xX a C T (^ ]a ) . (C.8) 
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